
namespace input
{
    std::string output_folder = "./res/data/";
    //! stepping function----------------------
    const int  step_type = 3;               // e指数增长（步长e指数增长）
    const double THz_max = 120;              //THz
    //! temporal-----------------------------
    const int    NC_pulse1   = 16;    //number Of Cycles
    const int    SP_pulse1   = 64;    //number of sampling points in one cycle
    const double Time_0  = 0.;
    const double phase0 = 0.;
    
    //! space_radial---------------------------
    const double Radius_0    = 0.;
    const double Delta_r_min = 5.0;   // nm    
    const double Delta_r_max = 200.0; // nm  

    //! space_z---------------------------------
    const double z_0      = 0.;    // Propagation direction Initial position
    const double distance = 4.e5;  // nm
    const double z_factor = 16;

    //! input laser1----------------------
    const double wavelength1 = 800.;  // nm 
    const double laser_intensity = 1.e14; // W/cm^2 
    const double FWHM   = 10.;     // fs
    const double phase1 = 0.;      //CEP
    const double waist1 = 5.e3;    // nm  
    const double focus1 = 2.e5;    // nm
   
    //! input laser2-------------------------
    const double E_factor = 0;   // E2/E1;倍频和基频场强比
    const double wavelength2 = 0.5 * wavelength1;  //nm 
    const int    NC_pulse2    = 2 * NC_pulse1;     //number Of Cycles
    const double phase2 = 0.;
    const double waist2 = 5.e3;  // nm 
    const double focus2 = 2.e5;  // nm
    
    //!medium------------------------------
    //type = frespace
    //index = air (1.0)
    //pressure = 1
    const double eta2 = 5.6e-19;                //Kerr;1.0e-19cm^2/W 
    const double density_of_neutrals = 2.7e19;  // 2.7e19cm^-3  5.4e18
} 



